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DETERMINATION OF LONGITUDINAL STABILITY AND CONTROL CHARACTERISTICS 
FROM FREE-FLIGHT MODEL TESTS WITH RESULTS AT TRANSONIC SPEEDS 
FOR THREE AIRPLANE CONFIGURATIONS 


By Clarence L. Gillis and Jesse L. Mitchell 


SUMMARY 

A test technique and data analysis method has been developed 
for determining the longitudinal aerodynamic characteristics 
from free-flight tests of rocket-propelled models. The technique 
makes use of accelerometers and an angle-of -attack indicator to 
permit instantaneous measurements of lift , drag , and pitching 
moments. The data , obtained during transient oscillations 
resulting from control-surface disturbances , are analyzed by 
essentially nonlinear direct methods (such as cross plots of the 
variation of lift coefficient with angle of attack) and by linear 
indirect methods by using the equations of motion for a transient 
oscillation. 

The analysis procedure has been set forth in some detail and 
the feasibility of the method has been demonstrated by data 
measured through the transonic speed range on several airplane 
configurations. Where comparisons were possible , values of 
aerodynamic parameters determined by different methods showed 
generally very good agreement. In addition , it was shown that 
the flight conditions and dynamic similitude factors for the tests 
described were reasonably close to typical full-scale airplane 
conditions. 

INTRODUCTION 

Aerodynamic testing by means of free-flight techniques 
offers several obvious advantages; thus, considerable develop- 
ment of suitable test and analysis methods has been made by 
many organizations. Free-flight techniques embrace a wide 
range of vehicle sizes from small gun-launched projectiles to 
full-scale airplanes. The technique described herein was 
developed for the purpose of permitting investigations of the 
longitudinal stability, control, and performance character- 
istics of aircraft by means of rocket-propelled models of 
fairly large size. 

Reference 1 contains a partial description of the technique 
for longitudinal stability and control measurements on free- 
flight models. The purpose of the present report which is 
based on references 2 to 4 is to describe the test procedure 
and to present a more comprehensive discussion of the 


method of analysis as it has been developed from a large 
number of model flight tests. For illustrative purposes, data 
are presented for three airplane configurations at transonic 
speeds. The model flight tests were performed at the Langley 
Pilotless Aircraft Research Station at Wallops Island, Va. 

SYMBOLS 

The axes system used is shown in figure 1. 
a n normal accelerometer reading, g units 
a t longitudinal accelerometer reading, g units 
a t transverse accelerometer reading, g units 
A aspect ratio, b 2 /S 
a logarithmic damping constant 

b wing span, ft 

Cyi cycles to damp to one half amplitude 
C c chord force coefficient 



Figure 1 . — System of axes. Each view presents a plane of the axes 
system viewed along the positive direction of the third. Angular 
displacements as shown are positive. 
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drag coefficient 
lift coefficient 
normal-force coefficient 
pitching-moment coefficient 
chord, ft 


mean aerodynamic chord 


CL 


c 2 dy , ft 


x 

Va 

y 

y' f 

yep 


section lift coefficient 
drag, lb 

length of cylindrical section of fuselage, in. 
acceleration due to gravity, 32.2 ft/sec 2 
moment of inertia about F-axis, slug-ft 2 
Iy 
qSc 

drag-due-to-lift constant 
constant used in aeroelasticity analysis 
radius of gyration about F-axis, ft 
lift, lb 

longitudinal distance between two normal accelerom- 
eters, ft 

tail length (distance from center of gravity to aero- 
dynamic center of tail), ft 
Mach number 
moment about A-axis, ft-lb 
moment about F-axis, ft-lb 
moment about Z-axis, ft-lb 
mass, slugs 
mV 
qS 

normal force, lb 

period of oscillation, sec 

calibrated static pressure, Ib/sq ft 

free-stream static pressure, lb/sq ft 

standard sea-level static pressure, 2116 lb/sq ft 

dynamic pressure, lb/sq ft 
Reynolds number 
radius, in. 
wing area, sq ft 

time to damp to one-half amplitude, sec 
time, sec 
velocity, ft/sec 
weight, lb 

displacement of accelerometer from center of gravity 
along X-axis, ft 

displacement of air-flow indicator vane from center of 
gravity along X-axis, ft 

distance along X-axis from any reference point, ft 
displacement of accelerometer from center of gravity 
along F-axis, ft 

spanwise distance along F-axis from plane of symmetry, 
ft 

spanwise distance along } -axis from fuselage side, ft 
spanwise position of center of pressure of exposed 
semispan, ft 


2 a displacement of accelerometer from center of gravity 
along Z-axis, ft 

a angle of attack, deg or radians 
«lo angle of attack at zero lift 

/3 angle of sideslip, deg 

7 flight-path angle, radians 

8 stabilizer deflection, deg or radians 

e downwash angle, radians 

„_j/L 

^ 6/2 

0 angle of pitch in equations of motion, radians; wing 

twist angle in plane parallel to plane of symmetry in 
aeroelastic analysis, radians 
A sweepback angle, deg 

jj. relative density factor, m/pSc 

p atmospheric density, slugs/cu ft 

^ roll angle, radians 

yaw angle, radians 
12 phase angle, radians 

co frequency of pitch oscillation, sec -1 

Subscripts: 

eg at center of gravity 

nose at model nose 

1 indicated 

ac aerodynamic center 

t trim, when used with a and C L 
T total 

r rigid 

e elastic 

0 at zero a and 8 

min minimum 

a accelerometer location 

aero aerodynamic 

b wing normal-force balance 

w exposed wing 

av average 

A single dot over a symbol indicates the first derivative 
with respect to time; a double dot indicates the second deriva- 
tive with respect to time. The symbols a, 8, and C L 
used as subscripts indicate the partial derivative of the 
quantity with respect to the subscript, such as 


C L 


a ‘ 


dCo 

da 


whereas the symbols a and q used as subscripts represent 
the partial derivative of the quantity with respect to the 

~c Oc 

uondimensional quantities — y and respectively, such as 



TEST PROCEDURE 

The rocket-propelled models are launched from the ground 
and are propelled by solid-fuel rockets. Information on the 
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steady-state and transient motions of the models is obtained 
by applying intermittent disturbances, by means of a con- 
trol-surface deflection, for example, and recording the re- 
sulting motions. These motions are measured by internal 
instrumentation in the models and are transmitted by tele- 
meter to ground receiving stations. Flight-path informa- 
tion is obtained by means of ground-based radar techniques. 
Data reduction and analysis of the telemeter and radar 
records yields data on the lift, drag, and pitching-moment 
characteristics of the models. Continuously recording 
equipment is used to permit the detection of rapid changes 
in the measured quantities. 

MODELS 

The four models used for illustrative purposes in the 
present investigation were part of a general research study 
of the effects of wing plan form on the transonic aerodynamic 
characteristics of airplane configurations. Sketches of the 
models are shown in figure 2 and photographs, in figure 3. 
The models included tfiree wing plan forms: unswept, 45° 
swept back, and 60° delta. The geometric characteristics of 
the various wings are given in table I. One model of each 
configuration was constructed with a solid-steel wing and 
an additional model of the swept-wing configuration was con- 
structed of solid duralumin to permit measurement of 
aeroelastic effects. 

The fuselages were bodies of revolution having cylindrical 
center sections. The nose and tail sections were defined by 
the ordinates of table II. The fuselages were composed of 
five separable sections constructed of sheetmetal skin at- 
tached to connecting rings. The vertical-tail surfaces were 
of constant thickness and had wedge leading and trailing 
edges. The vertical tails were constructed of wood and 
plastic and had aluminum surface inlays. 

The solid-aluminum horizontal tails were the same for all 
models and had plan forms and airfoil sections identical to 
those of the unswept wing. The horizontal tails were 
mounted on ball bearings built into the vertical tails (fig. 4) 
and were pivoted about a hinge line at 42 percent of the 
tail mean aerodynamic chord. For the swept-wing models 
the gaps at the roots of the horizontal tails were sealed by 
means of wiper-type seals. For the other two models these 
gaps were left unsealed. In order to produce the longitudinal 
motions of the models during flight, the horizontal-tail sur- 
faces were operated by control-actuating mechanisms in the 
rear sections of the fuselages. For the unswept-wing model 
this mechanism consisted of a cam operated by an electric 
motor. In the other three models this mechanism consisted 
of hydraulic power systems programed by means of electric 
motors. The control systems produced approximate square- 
wave motions of the horizontal tails. In both systems the 
time intervals of the control motions were preset by adjust- 
ing the speed of the electric motors. 

The swept-wing and delta-wing models were equipped 
with an extra vertical tail below the fuselage as indicated in 
figure 2. This tail was geometrically identical to the upper I 




Figure 2.— General arrangement of models. All dimensions in inches. 


vertical tail and was included on the models to minimize 
the possibility of lateral motions which had been encountered 
on other model flights. 

INSTRUMENTATION 

All models contained NACA telemeters which provided 
continuous records of the measurements of the internal in- 
strumentation. The telemeter was located in the nose of 
the model and the transmitting antenna consisted of a wire 
imbedded in the leading edge of the vertical tail. Measure- 
ments made on all the models included normal and longitu- 
dinal acceleration near the model center of gravity, angle of 
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Figure 3. — Photographs of the three configurations. 


TABLE I 

WING GEOMETRY 



Unswept 

Swept 

Delta 

Aspect ratio - — 

3. 00 

4. 00 

2. 31 

Taper ratio 

0. 40 

0. 60 

0 

Area (including fuselage), sq ft 

2. 68 

2. 68 

3. 15 

Airfoil section (free-stream) __ 

Hexagonal 

XACA 

NACA 



65A006 

65A003 

Thickness _ _ - - - 

0. 045c 

0. 06c 

0. 03c 

Mean aerodynamic chord, in._ 

12. 03 

10. 00 

18. 70 


attack, horizontal-tail deflection, total pressure, and a cali- 
brated static pressure. On the swept- and delta-wing models, 
measurements were made of transverse acceleration and wing 
normal force. In addition, the measurements for the delta- 
wing and flexible swept-wing models included a measure- 
ment of normal acceleration in the nose, which in conjunction 
with the normal acceleration at the center of gravity 
provided a measure of pitching angular acceleration. A gen- 
eral view of the telemeter installation and instrumentation 


TABLE II 

FUSELAGE NOSE AND TAIL ORDINATES 



x , in. 

r, in. 

0 

0. 168 

0. 060 

. 182 

. 122 

. 210 

. 245 

. 224 

. 480 

. 294 

. 735 

. 350 

1. 225 

. 462 

2. 000 

. 639 

2. 450 

. 735 

4. 800 

1. 245 

7. 350 

1. 721 

8. 000 

1. 849 

9. 800 

2. 155 

12. 250 

2. 505 

13. 125 

2. 608 

14. 375 

2. 747 

14. 700 

2. 785 

17. 150 

3. 010 

19. 600 

3. 220 

22. 050 

3. 385 

24. 500 

3. 500 


*d= 40.00 in. for unswept-wing model; d = 42.00 in. for all other models. 

in a typical model is shown in figure 5 and a more detailed 
view inside the model center section is shown in figure 6. 

Angle of attack was measured by a free-floating vane-type 
instrument located ahead of the model nose. I his instru- 
ment is more fully described in reference 5. I he total range 
of this instrument is approximately ±15°. In order to per- 
mit measurements to somewhat higher positive angles of 
attack, the sting to which xhe vane was attached was mounted 
at a negative angle, with respect to the fuselage center line, 
of 6° for the delta-wing model and 5° for the flexible swept- 
wing model. 

The total pressure was measured by a tube located under 
the fuselage (figs. 2 and 3) and the static pressure was meas- 
ured at an orifice located on the top center line of the body 
at a station about 0.6 body diameters back of the beginning 
of the cylindrical section. Previous to the flight tests of the 
models described herein, a special instrumentation test model 
was flown to calibrate these pressure-measuring locations. 
This calibration model (fig. 7) was composed of the forward 
portion of the fuselage of the airplane models with an angle- 
of-attack vane, total- and static-pressure orifices, and a longi- 
tudinal accelerometer, mounted on a two-stage rocket-pro- 
pelled test vehicle. Results from this test indicated that, 
within the accuracy of the measurements, the underslung 
total-pressure tube on this fuselage measured free-stream total 
pressure at small angles of attack up to at least 2.1. 
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(a) 



(a) Unsealed gap. 
(b) Sealed gap. 


Figure 4. — Empennage arrangement. 

The calibration of the static-pressure orifice from the in- 
strumentation test model is shown in figure 8. The static- 
pressure position error at supersonic speeds is small and rela- 
tively constant. The most important feature of the calibra- 
tion curve in figure 8 is the sharp change in the pressure at a 
Mach number between 0.97 and 0.98. Many subsequent 
checks of this pressure variation on other models having the 
same fuselage indicate that, within the accuracy of measure- 


ment, the sharp break always occurs at the same Mach num- 
ber. Thus, when this pressure break appears on the flight 
time history, the Mach number can be established at that 
point irrespective of the absolute accuracy of the pressure 
measurements. 

Wing normal force was measured by mounting the wing 
on a beam-type balance. (See fig. 9.) The balance was de- 
signed to permit vertical translation of the wing with respect 
to the bulkhead without angular motion. This vertical 
motion was measured by an inductance gage calibrated to 
give measurements of wing normal force. On the delta-wing 
model a strain gage was also attached to the wing-balance to 
furnish a duplicate measurement. The wing and balance 
assembly were fastened to the fuselage at the bulkhead. 
The gaps around the wing root where it entered the fuselage 
were sealed with rubber tubing or fabric. 

Ground-based instrumentation was also utilized during the 
model flights. A CW Doppler radar unit was used to meas- 
ure model velocity and a modified SCR 584 radar unit was 
used to obtain model position in space. An approximate 
measure of model rolling velocity was obtained from spin- 
sonde radio receiving equipment (ref. 6) which recorded the 
rolling rate of the telemeter antenna radiation pattern. 
Atmospheric conditions for the model flights were obtained 
from standard radiosondes released immediately after each 
flight. Fixed and manually operated 16-millimeter motion- 
picture cameras were used to record the launching and initial 
portion of each flight. 

PREFLIGHT MEASUREMENTS 

The model flights described herein took place at low alti- 
tudes, the portions useful for data analysis occurring at less 
than 10,000 feet. The resulting high dynamic pressures in- 
troduced the possibility of aeroelastic effects. Estimations 
based on static loading tests and calculations indicated that 
for the unswept and delta wings the aerodynamic effects of 
elasticity were negligible. For the swept wings, however, 
the effects were appreciable. Shown in figure 10 are the 
structural influence coefficients determined from static load- 
ing tests of the duralumin swept wing. The coefficients are 
shown as the twist in the free-stream direction per unit load 
applied at various spanwise stations along the 25- and 50- 
percent chord lines. 

For use in analyzing any flutter or buffeting vibrations 
that occurred in flight, the vibration frequencies of the 
model components were determined by suspending the 
models in shock cords and vibrating by means of an electro- 
magnetic shaker and by striking the wings and fuselages. 
The resonant frequencies were determined by visual obser- 
vation and from the records obtained from the model in- 
strumentation and telemeter during the vibration tests. 

After installation of all instrumentation, the weight, 
center of gravity, and moments of inertia of each model 
were measured. These quantities are given in table III. 
Moments of inertia were determined by swinging the models 
as pendulums and recording the frequency of oscillation. 
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Figure 7. — Test model for calibration of airspeed system, 


FLIGHT TESTS 

The models were launched from zero-length launchers; a 
typical launcher is shown in figure 1 1 . Each model was 
accelerated to maximum velocity by a single 6-inch-diameter 
solid-propellant ABL Deacon rocket motor. The rocket was 
fitted with an adaptor on the forward end into which the 
rear end of the model fitted and with a set of cruciform fins 
at the rear end to provide aerodynamic stability for the 
model and booster combination. Small fins, with a fixed 
angle of incidence, were mounted on the forward ends of 
some of the rockets. At separation of the booster from the 
model these fins cause the booster to be deflected away from 
the flight path of the model and this deflection aided in the 
radar tracking of the model. 

After rocket motor burnout the models separated from 
the boosters because of their lower drag-weight ratios com- 
pared with those of the model and booster combinations. 
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Figure 8. — Variation of the static-pressure ratio with Mach number 
for the fuselage orifice. 



Figure 9. — Wing balance components. 
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(a) Load applied along the 50-percent chord line. 

(b) Load applied along the 25-percent chord line. 

Figure 10. — Twist in the free-stream direction per unit load applied 
at various stations along the span of the duralumin wing. 


TABLE III 

MASS CHARACTERISTICS 


Unswept 

wing 

Swept 

Delta wing 

Steel 

wing 

Duralumin 

wing 

Weight, lb 

126 

146 

119 

142 

Center of gravity, 





percent c 

... 12.4 

24. 8 

4. 6 

20. 6 

/r, slug-ft 2 

___ 8.91 

10. 30 

10. 95 

11. 56 


The models contained no internal rocket motors and thus 
decelerated continuously after separation. The data used 
for analysis were obtained during the decelerating flight. 
The measurements made by the model internal instrumenta- 
tion were transmitted by the telemeter to ground receiving 
stations. A portion of the telemeter record for the triangular- 
wing model, showing all quantities recorded, is given in 
figure 12. Disturbances in the form of abrupt stabilizer de- 
flections between two limit positions were applied by the 
model control systems at predetermined time intervals during 
flight. The model motions during the short-period oscilla- 


tions following these disturbances were analyzed to deter- 
mine the aerodynamic characteristics. 

DATA REDUCTION AND CORRECTIONS 

The ground-based instrumentation is used to obtain a 
record of model velocity, Mach number, and dynamic pres- 
sure as a function of time during the model flights. The GW 
Doppler radar unit measures the velocity component of the 
model in a radial direction from the radar location. The 
models are launched from a point as close as possible to the 
radar unit so that the velocity so measured is essentially 
the total velocity. Model-flight-path information in the 
form of elevation, azimuth, and range as a function of time 
is obtained from the SCR 584 radar unit and is used to cor- 
rect the velocity obtained by the CW Doppler radar unit to 
total velocity. Atmospheric conditions obtained from 
radiosondes are then used with the corrected velocity to 
obtain Mach number and dynamic pressure. 

For the models discussed herein no data from the CW 
Doppler radar unit were obtained after separation of the 
models from the boosters, and for the unswept wing model 
no data from the SCR 584 radar unit were obtained after 
this point. Failure to obtain such data was due to the rapid 
maneuvers of the models which made tracking difficult. 

When no velocity data from radar units are available, 
the Mach number is obtained from the standard compress- 
ible-flow equations utilizing telemetered total pressure and 
a static pressure obtained from radar-measured altitude and 
radiosonde pressures. When no radar-measured altitude 
information is available, the static pressure is obtained from 
the calibrated static orifice on the model as described 
previously. 

The oscillograph traces of the telemetered measurements 
are reduced to physical quantities by use of the calibrations 
of the recording equipment made immediately after each 
(light and the calibrations of the individual instruments 
made as near as possible to the time the instruments were 
installed in the model. The measured accelerations are 
then reduced to coefficient form for data-analysis purposes. 

Corrections are usually required to some of the measured 
quantities to obtain the desired information. The acceler- 
ometers used to measure accelerations at the center of gravity 
can not all be located on the center of gravity, and corrections 
may be required to account for the distance by which the 
accelerometers are displaced from the desired point. These 
corrections must account for the effects of angular velocities 
and angular accelerations. For accelerometers located at 
distances x a , y a , and z a from the center of gravity (distances 
positive in the direction of positive A"-, F-, and Z-axes of 
fig. 1), the equations for correcting accelerations to the center 
of gravity are: 

o*,e*=a».<+“ [— x tt (ij— 'Pvd+yaiv+bi')— z a (^ 2 +0 2 )] (!) 
«/.<*=«/,( + “ [— Xait + Op) + l/airpt+V*) +Z a (v — ' fy)] (2) 
a,, «f=Oi.<+- [ar a (tf 2 +ie)+2/a(^— z«(0+^p)] (3) 
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(a) Before launching. 


Figure 11. — Model launching. 


(b) After launching. 



Wing load Body static pressure Normal acceleration at center of gravity Free -stream total pressure 



Normal acceleration in nose 




Longitudinal acceleration 

\ 

Horizontal-tail position 



1 - r r 


Wing load Transverse acceleration Model angle of attack 



Figure 12. — Portion of telemeter record for delta-wing model. Mach number approximately 1.1. 
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It is obvious that, in general, the corrections to the 
accelerometers are very complicated. The desirability of 
eliminating lateral motions is indicated since the corrections 
are thereby considerably simplified as follows: 


'n, eg ^ n , (d “ ( X a 6 ) 

(4) 

& l,cg (*^a^ 2 Za.0 ) 

(5) 


The pitching velocity 6 and angular acceleration 6 may be 
approximated in the manner described subsequently in 
conjunction with analysis of the pitching-moment data. 

The air-flow direction vane on the nose of the model 
measures angle of attack plus an increment in air-flow 
angularity caused by pitching velocity. The corrected angle 
of attack (for a and a f in degrees) is given by 

a:=a l + 57.3 pr 6 (6) 

The wing normal forces obtained from the wing balance 
include not only aerodynamic forces but also inertial forces 
exerted by the wing and the moving parts of the balance. 
The total normal forces read by the balance are corrected for 
these inertial forces to yield the aerodynamic normal force 
as follows: 

N a ero = N T -\-a„W\ (?) 

The resulting aerodynamic normal force is converted to 
coefficient form and the exposed wing lift coefficient is 
assumed to be given by C LtW =C NtW cos a. 

Because the frequency-response characteristics of the 
instruments used in these tests were satisfactory, no correc- 
tions were necessary for instrument lag or amplitude attenu- 
ation effects. 

METHOD OF ANALYSIS 

Many mathematical procedures have been developed for 
analyzing transient oscillation data (refs. 7 and 8) to obtain 
aerodynamic derivatives. The method of analysis described 
herein makes use of measurements obtained during the 
transient oscillations resulting from step-function disturb- 
ances. The procedure adopted as a practical method of 
analysis was to obtain as much information as possible from 
direct measurements and to avoid extensive mathematical 
manipulation of the data. Thus problems of handling such 
questions as nonlinearities and the form of the derivatives 
could be minimized, and exploration of regions such as the 
stall could be conducted with little restriction. Successful 
application of the procedure depended upon the development 
of an angle-of-attack indicator with the required accuracy 
and frequency-response characteristics. This was success- 
fully accomplished and is described in reference 5. Determi- 
nation of some of the aerodynamic characteristics required 
a solution of the equations of motion, however, and this 
solution as well as the use of direct measurements will 
be described. 

DIRECT MEASUREMENTS 

Lift and drag. — The accelerometer measurements at the 


center of gravity provide normal-force and chord-force 
coefficients according to the following relations: 


„ W 
Cat a„ qS 

(8) 

and 


n w 

(9) 

The total lift and drag coefficients are then 

obtained from 

C L =C N cos a—C c sin a 

(10) 

C D =C c cos a-\-C N sin a 

(ID 


The lift characteristics are determined by plotting lift 
coefficient against angle of attack. The first one or two cycles 
of a transient oscillation are generally used for this purpose 
to obtain the largest amplitudes. Such a plot may also be 
made during a nonoscillatory motion. The change in Mach 
number during one cycle of an oscillation on the flight tests 
reported herein was about 0.02 at supersonic speeds and 
about 0.007 at subsonic speeds. Thus, unless the aerody- 
namic parameters vary very rapidly with Mach number, the 
error involved in assuming that the Mach number is constant 
should be small. Figure 13 presents typical plots of lift 
coefficient as a function of angle of attack. These data are 
for the 45° swept-wing configuration (duralumin wing model). 
With the exception of the results at a Mach number of 0.92, 
these data were obtained for the angle of attack both increas- 
ing and decreasing with time. For this particular model the 
effect of rate of change of angle of attack with time is not 
detectable from the data. At a Mach number of 0.92, data 
on the variation of lift coefficient with angle of attack were 
obtained at angles of attack up to 18° when the model per- 
formed a pitch-up maneuver for the horizontal-tail deflection 
of —4.60°. From results such as these, the lift-curve slope 
(including any nonlinearities) as well as angle of zero lift, 
maximum lift coefficient, and the lift phenomena associated 
with the separated flow at high angles of attack may be 
determined. 

The lift coefficient as determined is, of course, the total lift 
coefficient and includes any contributions from derivatives 
such as C Lq and C L& . Normally, the effects of these con- 
tributions are negligible as indicated by the data of figure 13 
but they may be calculated by using estimated values of 
C Lq and C Lit . Such calculations have been made for several 
cases and it has generally been found that, because of the 
phase relations between a and q, the effect of C Lq and C Lit 
is to cause a slight displacement of the curves of (\ against 
a in such a way as to indicate a phase lead of the lift coeffi- 
cient compared with the angle of attack but with no measur- 
able change in the lift-curve slope. 

The data from the wing normal-force balance are also 
analyzed by the direct method. Results for the swept-wing 
configuration are typical of the data obtained. An examina- 
tion of the results (fig. 14) indicates that these data also fur- 
nish the variation of lift coefficient with both angle of attack 
and Mach number. 
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Figure 13. — Lift data for duralumin swept-wing model. 
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The drag information is generally of most use when plotted 
in the form of polars of C D against Cl . Typical results are 
illustrated in figure 15 for the swept-wing configuration. 
From such plots the minimum drag coefficient and the varia- 
tion of lift-drag ratio with lift coefficient may be determined. 
Since two control deflections are available, some measure of 
the increment in drag due to trimming the airplane can also 
be obtained. Note that only for the drag curve at a Mach 
number of 0.97 is the effect of the small change in Mach 
number significant. The oscillation represented by these data 
occurred during the rapid transonic drag rise. 

The variation of drag with lift is sometimes represented by 
an equation of the form 

C D =C D ,mln+KC L 2 ( 12 ) 


where A plot of C D against C L 2 will yield a linear 

curve whose slope is K if the preceding equation adequately 
represents the data. For a wing with camber or for any con- 
figuration change which produces an effective camber, 
the lift coefficient for minimum drag will generally not be 
zero, and such a drag curve may frequently be represented 
by the parabolic equation 

C D =C D , min +K(C L -C L ,c D ^ in y (13) 


for the purpose of determining K. 

Pitching moment.— The total pitching-moment coefficient 
is obtained from the relation 


c m = 



(14) 



Figure 15. — Drag data for duralumin swept-wing model. 



Figure 15. — Concluded. 


The pitching acceleration may be measured by a pitching 
angular accelerometer or, as was done for some of the model 
results presented herein, by two normal accelerometers dis- 
placed along the longitudinal axes. With one accelerometer 
in the model nose and one at the center of gravity, for 
example, the pitching acceleration will be given by 


e=2^ 



(15) 


This result can be easily obtained from equation (4) if the 
2 a -distance of both accelerometers is the same. If measure- 
ments of 0 are not available, an approximation can be found 


from the relations 



(16) 


6=a-\- 7 


(17) 


The quantity a is obtained by differentiating the measured 
a-curve and y is calculated from the equation 


y=y (fln.cg— COS e COS <f>) 


(18) 


DETERMINATION OF LONGITUDINAL STABILITY AND CONTROL FROM FREE-FLIGHT MODEL TESTS 


13 


For the flight conditions existing on the rocket-propelled 
models (large q and T 7 ), the last term in the latter equation 
is very small and may generally be neglected; thus the need 
for measuring 6 and p is obviated. 

The total pitching-moment coefficients obtained as de- 
scribed above contain contributions due to the damping 
terms C mq and C ma and, unlike the corresponding lift terms, 
these are not usually negligible although they will generally 
be small compared with the pitching moments due to a . 
The pitching moments due to the damping terms may be 
subtracted from the total pitching moments to obtain 
pitching moments as a function of angle of attack by 

C m {a)=C*~ C mq y (19) 

where a and y are obtained as described previously. The 
sum of the damping derivatives C mq -\-C m6i for use in this 
calculation is obtained as described later. The individual 
value of C mq cannot be determined from the data but must 
be estimated for the purposes of applying this correction. 

The resulting pitching-moment coefficients may be plotted 
against angle of attack or lift coefficient to obtain the static 
stability derivatives C ma and dCJdC L , including nonlinear 
effects, the trim values of a and C L for a given control de- 
flection, and the value of C m , o- From the increments in 
C m for different control deflections at a constant angle of 
attack, a measure of the control effectiveness C mh may be 
obtained. Typical pitching-moment curves are shown in 


figure 16 for the swept-wing configuration. These pitching- 
moment data were obtained from two linear accelerometers 
and the results have been corrected for the previously dis- 
cussed pitching moments due to the damping terms. A 
significant effect of Mach number during one oscillation is 
evident only for the curve at an average Mach number of 
0.97. In this region there is apparently a rapid change in 
pitching moment with Mach number at a constant lift 
coefficient. The data at a Mach number of 0.92 are typical 
of results obtained when the configuration has a highly non- 
linear variation of pitching-moment coefficient with lift 
coefficient. Note that even relatively smaller nonlinear 
variations, for instance at M= 0.94, are readily apparent 
when the data are presented in this manner. 

It is possible to analyze the data obtained from two 
accelerometers without first reducing the measurements to 
coefficient form by combining equations (14) and (15) to 
obtain an equation for the static stability slope dC m ldC L as 
follows: 

< ^k=~ hr. (l (20) 

dC L TflCl \ (ML ntC g / 

Although the slope so determined is the rate of change of the 
total pitching-moment coefficient, it should be essentially 
unaffected by the damping moment if the slope is measured 
at C m = 0 because at this point the rate of change of the 
damping moment is near a minimum. 

Control effectiveness. — If the amplitude of the oscillatory 
motions is small at the time the control-surface deflection is 



(a) Horizontal-tail deflection -4.60°. (b) Horizontal-tail deflection -0.72°. 

Figure 16. — Pitching-moment data for the duralumin swept-wing model as obtained from the two-accelerometer method. Flagged 

symbols indicate decreasing Cl. 
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initiated and the control movement is sufficiently rapid, the 
changes in lift and pitching moment due to control deflection 
can be detected and a direct measure of the control effective- 
ness parameters C mi and C L& can thus be obtained. For the 
models of the present investigation the motions were, in 
general, not of sufficiently small amplitude when the control 
motions were initiated to permit reliable values of these 
parameters to be determined by this method. 

TRANSIENT OSCILLATION ANALYSIS 

The method of analysis used herein applies to the free 
oscillation resulting from a disturbance. For the models 
considered herein, the disturbance consisted of a rapid 
control-surface deflection; however, other types of disturb- 
ances may be used, such as the impulse from a small rocket. 
The complete derivation of the equations used will not be 
given herein because it is fairly simple and may be found in 
a number of sources. Only the final results and the method 
of applying them to free-flight models will be shown. The 
discussion has been kept general in character for purposes of 
application to other models. 

In order to simplify the analysis and to permit the deter- 
mination of equations for the more important aerodynamic 
derivatives, several assumptions are necessary. It is 
assumed that, during the time interval over which each 
calculation is made, the following conditions hold: the 
forward velocity and Mach number are constant; the aero- 

(Ict 

dynamic forces and moments vary linearly with «, 6 , 5, 
dB 

and -tt ; and the model is in level flight before the disturbance 
at 

is applied. A discussion of the effect of these assumptions 
upon the results is presented subsequently. 

The first assumption mentioned in the previous paragraph 
effectively limits the longitudinally disturbed motion of the 
aircraft to two degrees of freedom: translation normal to 

the flight path and rotation in pitch about the center of 
gravity. The equations of motion resulting from these 
assumptions are 

m ' ^“^) = CL,0+CL a «+Ctj5 (21) 

r ^=c m . 0 +c m «+c m , ± f t +c m , ~ f t +c m ,s ( 22 ) 

where all angular quantities are in units of radians. When 
these equations are solved, the following equation for the 
free oscillation of the angle of attack is obtained: 

oL—Ca at cos (cot 12) ~\~oc t (23) 

where C is a constant. The symbol a t denotes the steady 
state or trim angle of attack which will exist after the oscilla- 
tion has damped out and represents the mean value about 
which the angle of attack oscillates. The first term repre- 
sents the oscillation about the trim angle. Figure 17 is a 
schematic plot showing a typical record of the angle-of-attack 
response following a step deflection of the aircraft control 
surface. 

The constants a, c o, and a t in equation (23) are independent 
of the initial conditions and the analysis consists essentially 



/ 


Figure 17. — Typical angle-of-attack oscillation as used in analysis. 

of finding the numerical values of these constants from the 
measured data and from these constants determining the 
stability characteristics of the configuration tested. The 
constants C and 12 depend upon the initial conditions and are 
not used in the analysis; thus, their numerical values need 
not be known for the type of analysis considered herein. 
From the envelope curves enclosing the oscillations, the 


damping constant a can be determined, 
to the notation in figure 17 

If reference is made 

Aa 2 =Aoi l e a(t2 ~ tl) 

(24) 

i Aa 2 

,o6 ' -a 

a t 2 / 1 

(25) 

The constant co defines the frequency or 
tion and is given b}^ 

period of the oscilla- 

N In 

CM 1^1 

II 

3 

(26) 


The constant a, is simply the value of a after the oscillation 
has damped to a steady value or is the value of a on the mean 
line of the oscillation as shown in figure 17. 

In order to determine the constants a, co, and a t from the 
measured data, it is necessary first to fair envelope curves for 
the oscillation which should be logarithmic curves according 
to equation (23). The mean line between the two envelope 
curves is drawn and values of Aa 2 , A«i, P, and a t can then be 
determined and a and co can be calculated. 

The basic data obtained from this analysis of the oscilla- 
tions are usually plotted as a function of Mach number. The 
results of such an analysis for the delta-wing configuration 
are shown in figures 18 to 20. The solid lines on the trim 
curves in figure 18 were obtained from the mean lines of the 
oscillations resulting from the rapid control deflections. The 
period of the pitch oscillation in this case was determined by 
measuring the time difference for trim crossings. The data 
for the two different control deflections (fig. 19) define two 
distinct curves. Although some scatter of the data points is 
evident, it can be seen that a variation of period with Mach 
number can frequently be detected within each oscillation. 
The variation of amplitude ratio with time plotted on semi- 
log paper so that straight-line fairings of the data can be 
utilized is illustrated in figure 20 (a). The variation with 
Mach number of the slope of these curves, given as time to 
damp to one-half amplitude, is shown in figure 20 (b). The 
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(a) Trim angle of attack. 

(b) Trim lift coefficient. 

Figure 18. — Trim lift and angle-of-attack characteristics of delta-wing 

model. 
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(a) Variation of amplitude ratio with time. 

(b) Time to damp to one-half amplitude. 

Figure 20. — Damping-in-pitch characteristics of delta-wing model. 
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Figure 19. — Period of pitch oscillation of delta-wing model. 


constant a can then be obtained from these results as 


—log e 2_— 0.693 
T'A 


(27) 


The success of this procedure as illustrated depends in part 
upon the oscillations being rather lightly damped so that 
several cycles are available during each oscillation to permit 
the fairing of envelope curves. If the oscillations are heavily 
damped, other methods of analysis will be necessary. 

The analytical solution of equations (21) and (22), which 
include those aerodynamic derivatives which previous experi- 
ence has indicated have an important influence on the motion, 
shows that the constants a and co are given by 



or 


2m' [ r ' £ “ ' m °+ Cm a ) 2 Qy) ] 


1 

1 

- ~ c 1 


r 

2V\ 

a2 

, I’m’) a 


(29) 

(30) 


The aerodynamic derivatives C Lq and C L& have been omitted 
from this analysis. The effect of these and other omitted 
terms upon the results is discussed subsequently. Equations 
(28) and (30) may be rearranged to give 




_ 4 rv( , c L \ 

c \ 2m' ) 


( 31 ) 


<?„,„= -/'(co 2 +a 2 )- 


7 - C n CJ La 
2V m' 


( 32 ) 


It is necessary therefore to know the value of C La and C m to 
use in equations (31) and (32) for calculating static stability 
and damping. The lift-curve slope C La for use in equation 
(31) is found by plotting C L against a as described previously. 
Numerical calculations have shown that the last term in 
equation (32) will probably always be very small compared 
with the first term (less than 1 percent) and may thus be 
omitted. Its effect may be estimated in any case. If this 
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term is omitted, the static stability parameter is then 


C m = -/'( co 2 +a 2 ) (33) 

From C Lay a , co, and the mass characteristics of the model, 
the damping factor C mq +C m6i and the static stability deriv- 
ative Cm a can be calculated by use of equations (31) and (33). 
The results for the delta -wing configuration shown in figures 

21 and 22 are typical. Since it is obvious from figures 21 and 

22 that these quantities have a nonlinear variation with 
angle of attack, the values obtained from the linear analysis 
must be considered as an average or equivalent linear value 
of the derivatives. The aerodynamic-center location as a 
fraction of the aerodynamic chord is then 


fgc %cg 


C C 



(34) 


where x ac and x cg are measured rearward from the leading 
edge of the mean aerodynamic chord. 

The average longitudinal control parameters C m ,o and C mb 
can be determined from the following analysis of the trim 
and stability data. The equation for the steady-state angle 
of attack a t obtained by solving equations (21) and (22) 
with 0 and a equal to zero is given by the following equation: 



very small and can be neglected, particularly when the air- 
plane relative density is large. Their magnitude should, 
of course, be checked. When these terms are ignored, 
equation (35) can be expressed as 


( m ,0 , da t ? 

a ‘ = ~cr~cz 5=a ‘^ + ^s s 


(36) 


The values of a t are determined for each oscillation and are 
then plotted against Mach number and curves faired through 
the data as illustrated in figure 18. The slope da t jdb is 
obtained from the increments between the curves, and 
oi t , 6=o is obtained by interpolation between the two a* -curves. 
From equation (36), it is now possible to calculate C m> 0 aud 
C 

1 ni 5' 

An application of the above procedure to the results for 
the delta-wing configuration gives the values of C m> 0 and 
C m8 shown by the solid-line curve in figure 23. The average 
value of C Wa from figure 21 was used with equation (35) to 
obtain the curves shown. The values of C m& and C m> 0 as 
obtained from the previously discussed method utilizing the 
pitching-moment curves is shown for comparison. The 
agreement between the two methods of determining C mb and 
C m ,o is considered good in view of the nonlinearity of the 
basic pitching-moment data for this configuration. 

Several procedures may now be used to determine C Lb , 
The value of C Lb can be obtained from the increment in C L 
at any given angle of attack between the lift curves plotted 
for two successive oscillations at different control deflections. 


The terms in equation (35) involving 



—y are generally 
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Figure 21. — Variation of longitudinal stability parameter Cm a with 
Mach number. Delta-wing model. 
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Figure 22. — Variation of rotary-damping coefficient with Mach 
number. Delta-wing model. 
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(a) Control effectiveness. 

(b) Pitching-moment coefficient at a = 5«=0. 

Figure 23. — Analysis of trim characteristics. 
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In equation form, 


,, (A C L ) 

Ll s =- 


a — Constant 

A5 


(37) 


In addition, C Lb can he calculated from the equation 


/o doit ■ dC 


L, t 


db db 


(38) 


where dC Ltt /db is found from the C Ltt and b curves in the same 
manner as dajdb . Another check on the value of C Lb can be 


obtained from 


Cll — 


hfc 


(39) 


where l t fc is the longitudinal distance from the center of 
gravity to the center of pressure of the lift caused by control- 
surface deflection. For a conventional airplane configura- 
tion with a horizontal tail, this distance can be estimated 
fairly accurately. 

The three methods given for determining C Lb are not 
all equally applicable to all configurations. For an airplane 
with a conventional tail surface to the rear of the wing, the 
value of C Lb is usually small, and the distance to the center 
of pressure of the lift due to control deflection can be esti- 
mated fairly accurately. Equation (39) should then give 
the most reliable results for C Ly Conversely, for a configura- 
tion having all-movable wings near the center of gravity 
C Lb is of the same order of magnitude as C La and the distance 
to the center of pressure of the lift due to control deflection 
cannot be accurately estimated, since this lift includes an 
increment caused by wing downwash on the tail. For this 
case equations (37) and (38) should give a more reliable 
result. 

Figure 24 presents a comparison of the results obtained on 
the unswept-wing configuration by using the three methods. 
For this particular configuration it is not readily apparent 
from the data which method gives the most consistent results 
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Figure 24. — Comparison of three methods of determining tail lift 
effectiveness. Unswept-wing model. 


but from the previous discussion the results from equation 
(39) are considered to be the most reliable. 

EFFECT OF ASSUMPTIONS IN ANALYSIS 

The method of analysis described and illustrated herein 
has been developed by use of a number of simplifying as- 
sumptions, as noted in the analysis section. It is instructive 
to examine some of these assumptions to determine their 
validity where possible. 

One of the assumptions was that of constant forward 
velocity. The variable forward velocity which actually 
existed during the model flights may be considered to have 
several effects: (1) the small perturbation velocity which is 
present when the third degree of freedom of the longitudinal 
motion is permitted; (2) the variation in the mean velocity, 
consisting of an essentially constant deceleration over a small 
time interval; (3) the effect of unsteady air flow on the aero- 
dynamic characteristics of the various components of the 
aircraft; and (4) the variation of the aerodynamic deriva- 
tives with Mach number. 

The effects of the small perturbation velocities may be 
considered by comparing the calculated motions from a 
two-degree-of-freedom system of equations with those from 
a three-degree-of-freedom system. This comparison was 
made by using the characteristics of the unswept wing model 
of this report. Adding the third degree of freedom caused 
a change in oscillation period and damping of less than 1 
percent which is small compared with the overall test and 
data-reduction accuracies. 

The effects of a variation in the mean velocity have been 
considered by several investigators (refs. 9 and 10), and it 
has been found that a first-order effect of acceleration con- 
sists of the addition of terms dependent on V /V 2 to the damp- 
ing and frequency constants of the motion. A numerical 
check, using the experimental data from this report, indi- 
cates that the deceleration existing during the time the 
experimental data were obtained could cause a change in 
damping of less than 1 percent and a very much smaller 
change in period. 

Reference 11 investigates analytically the effect of acceler- 
ated air flow on the pressure drag and lift-curve slope of thin 
wings in the transonic and supersonic speed ranges. The 
wing is considered as decelerating from supersonic speeds, and 
this condition actually exists for the rocket-propelled models. 
The largest effect occurs at a Mach number of 1.0 where for 
unaccelerated air flow the aerodynamic quantities are infinite 
and for decelerating air flow the aerodynamic quantities are 
finite. Numerical calculations in reference 11 indicate that, 
at M=1.02, for example, the pressure drag of a wedge air- 
foil and the lift-curve slope of a flat-plate airfoil of the size 
on the models discussed in this report would be decreased by 
16 percent if the airfoil were decelerating at a rate of 376#. 
The actual deceleration obtained on the models herein was 
about 4# at supersonic speeds. This effect may therefore 
be considered negligible. 

Since the model is decelerating, the Mach number will 
change during each oscillation and the aerodynamic deriva- 
tives will also change. The values obtained in the analysis 
are thus average values over a small Mach number interval. 
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The effect of the varying Mach number on the calculated 
results is minimized by determining the required information 
over the shortest possible time interval. The lift-curve 
slope can be determined during one-half cycle of an oscilla- 
tion, for instance. The static-stability derivative C ma can 
be obtained by measuring each half-period and multiplying 
by two, plotting the results against Mach number, and mak- 
ing the computations from a faired curve as shown in figure 
19. The damping time can also be measured several times 
during an oscillation and plotted against Mach number. 
Similarly, the values of C Lt o, Cl.i, and a t can be determined 
for each oscillation and plotted against Mach number; the 
quantities a L0 , dC^tldb^ dajdb, and a*, 5=0 can be obtained 
from faired curves through the points rather than from the 
increments between measured points at different Mach 
numbers. The change in Mach number during the flight 
tests for one cycle of an oscillation was about 0.02 at super- 
sonic speeds and about 0.007 at subsonic speeds. r I bus, 
unless the aerodynamic parameters vary very rapidly with 
Mach number, the error involved in assuming that the Mach 
number is constant should be small. 

The method of analysis is not strictly valid if the aero- 
dynamic derivatives are not constant at a given Mach 
number. The primary effects of such nonlinearities can be 
determined, however, by choosing control-surface deflections 
which cause the model to oscillate over different ranges of 
angle of attack. This effect was evident on two of the 
model flights described in this report where different values 
of C La and C ma were obtained for positive and negative con- 
trol deflections in the transonic region. If the aerodynamic 
derivatives are extremely nonlinear within the region covered 
by one oscillation, then the values obtained in the analysis 
are only average values indicative of the trend of the data. 

Certain of the derivatives can be nonlinear without 
seriously affecting the results. For example, the period of 
the oscillation is almost completely determined by C ma , all 
other terms having only a very small influence; thus, non- 
linearities in all other derivatives would not appreciably 
affect the calculated values of C ma . Considerable judgment 
is necessary in interpreting the data when evidences of non- 
linearities exist and other more laborious methods of analysis 
may be necessary. 

When the three-degree-of-freedom calculations were made, 
certain other effects were also investigated. These were the 
effects on period and damping of the omission of C Lq and 
C L . from the analysis and the effects of the assumption of 
zero mean flight-path angle. For the flight conditions of 
the models discussed herein, these effects were negligible. 

AEROELASTIC ANALYSIS 

The method used to estimate the change in lift and center 
of pressure due to elastic deflection of the swept wings is a 
form of strip theory. The information needed to apply 
this method is data on twist along the wing due to unit loads 
applied along the centers of pressure (assumed to be along 
the 0.25 and 0.50 chord line), an assumed rigid-wing span- 
loading curve, and an estimate of the rigid-wing lift-curve 
slope. 


The portion of the half wing which is elastic (that portion 
outside of the fuselage) is divided into as many parts as 
desired. The rigid lift per unit angle of attack of these 
sections is then estimated from the span loading. The 
increment in lift due to twist at the section is then assumed 
to be given by this lift-curve slope multiplied by the twist. 
Since the twist at a given section depends on the final equilib- 
rium lift at all sections, a set of n simultaneous equations 
(equal in number to the number of sections the wing is 
divided into) must be solved; this solution gives the resultant 
elastic lift at each section. The equations are set up by 
use of the ratio of the elastic lift at each section to the rigid 
lift at the section. The independent variable used in the 
equations is the rigid lift-curve slope times the dynamic 
pressure C L<xr q. 

Drawings for the exposed-wing panel and span-loading 
curve divided into five sections are given in sketches 1 and 2. 
The equations for each section of the exposed-wing panel 


are 


L X =K'Q ^ f 1 f~ ^ ~dr](oir + Am) 

^ J V. , \' L^ac / rigid 


fuselage 
side 


(40) 



Sketch 1. 
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where the Z-terms are the loads on each section of the 


exposed wing, 77 K' 


1 


Q is defined 


/2 "7 ( 7 ) * 

^fuselage ' 
side 

as C Lot r w q, and K lf K 2 , . . . K 5 are defined as follows: 


K^K' 


K : =K' 



(41) 


The ratio of the twist at each section to the rigid angle 
of attack can be derived from equations (40) and (41) for 
the following expressions: 


Aq?i L x 

a r a r K\Q 


(42) 


Aa 5 _ L 5 j 

a r a r K b Q 


r Phe influence coefficient B tj is defined as the twist at station 
i due to a unit load applied at station j, and the equations 
for the twist of each wing section due to loads L iy L 2 , . . . 
L 5 are as follows: 


variable: 


— (hi-jrn)+— *»+— «u+— «u+— 0..— 1 

a r \ Oir ol t a T Oir 

— 0 «+— 0 M +- $ u +— ®m+— jZj)—! 

CL r Oi T Oi T Oi r Oi T \ 

After equations (44) are solved for Zi/« r , L 2 \a T 
L h \a T for each value of Q, 


(44) 




( L\ - L 2 | Lz j L\ | ZA 
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for the exposed wing. 

The effect of wing-inertia loading acting in opposition to 
aerodynamic loading has been neglected in equations (44). 
The addition of a wing-inertia term for a case where the 
wine; is divided into two sections is illustrated in sketch 3: 



— =— 0U+— 012 +— 0,3+— * 0,4+— I 
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When equations (42) are substituted into equations (43) 
for A ajar, A a 2 /a r . . . A a 5 /a r , the following set of five 
simultaneous equations is obtained with Q as the independent 
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where 

W total model weight 

Wi weight of inboard section of exposed half-wing 
W 2 weight of outboard section of exposed half-wing 
L t total lift of model 
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The inboard movement in the spanwise center of pressure 
which results from the reduction in lift due to aeroelastic 
wing twist causes a change in aerodynamic-center position. 
The lateral center of pressure of the exposed elastic wing can 
be found from L\\a T , L 2 lot T , . . . L b /a r solved from the 
simultaneous equations (44) as follows: 


/ / x _Liy[+L 2 y2+L 2 y^+L A 7y A -\-L 5 7j' 5 

KVcv) e -~ L 1 +L 2 +L,+L a +L 5 


(47) 


Measurements of the lift-curve slope of the steel and dural- 
umin swept-wing models were extrapolated to rigid values of 
lift-curve slope by the method indicated previously. These 
values of rigid wing lift-curve slope and the measured 
influence coefficients of figure 10 were used with equations 
(44) to calculate the ratio of the elastic lift-curve slope to the 
rigid lift-curve slope. The agreement between t he calcu- 
lated and measured results shown in figure 26 indicates the 
validity of the aeroelastic analysis method described. 


If the span load is assumed to move inboard along the 
half- or quarter-chord line, the increment in aerodynamic- 
center position of the complete model due to the inboard 
movement can be found from the following expression: 


dC m _ T (2/cp) r — (Vcy) 

i>C L l e 


-Tan A 



(48) 


For the tan A use either tan A c/4 or tan A c/2 whichever is 
appropriate. 

An application of this aeroelastic analysis to the swept-wing 
configuration illustrates the usefulness of the results. If the 
aeroelastic effect is represented by a one-point solution, the 
following result will be obtained : 



This simplified result is useful in determining rigid values of 
lift-curve slope from measured elastic values. An applica- 
tion of this result is given in figure 25. The results shown are 
for the steel and duralumin swept-wing models and an addi- 
tional point from some results (ref. 12) for a configuration 
having the same wing but constructed of wood and dural- 
umin. Two different values of (0/L) ref were used, one for a 
load applied on the 0.50 chord and the other for a load 
applied on the 0.25 chord. Straight lines through either set of 
points have a common point of intersection at qS(d/L) ref = 0. 

This may be taken as the rigid value of 
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ACCURACY 

The accuracy of the results from any one test cannot be 
precisely defined. The following discussion of the probable 
order of magnitude of the accuracy was obtained from a 
consideration of the accumulated experience from repeated 
calibrations of the instruments and recording equipment and 
of flight test experience where duplicate models were flown 
or duplicate measurements made by different instrumenta- 
tion. 

The absolute magnitude of a telemetered quantity will in 
most cases be accurate to within 1 percent of the total cali- 
brated range and the errors will very seldom exceed 2 percent . 
The incremental values or rates of change of the telemetered 
quantities should be accurate to 1 percent or better in almost 
all cases. For maneuvering models, such as those described 
herein, the radar equipment w r ill generally give results 
accurate to about 1 percent in the early portions of flight with 
errors greater than 1 percent in later portions. 

Based on the foregoing and other considerations, the 
estimated probable accuracy of some of the basic quantities 
involved in the tests is given in table IV. The accuracy 
of at for the unswept-wing model is considerably less than 
that for the other models because the accelerometer was 
calibrated to cover a much greater range of acceleration in 
order to include the accelerating portion of flight during 
boost as w'ell as the decelerating portion, whereas for the 
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Figure 25. — Variation of reciprocal of model lift-curve slope with wing 
stiffness at M«= 1.11. 


Figure 26. — Comparison between experimental and calculated elastic 
to rigid lift-curve-slope ratios for the duralumin swept wing. 
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other models the longitudinal accelerometers were calibrated 
to include the decelerating portion plus 1 g positive accelera- 
tion and were thus against stops during boost. The effects 
of the estimated accuracies of these basic quantities on some 
of the derived aerodynamic parameters are given in table V. 
In these tables certain quantities are included as ratios, such 


TABLE IV 

ESTIMATED ACCURACIES OF BASIC QUANTITIES 
[All increments may be positive or negative] 


Quantity 


Estimated accuracy 
for — 

M= 1.35 

Af=0.9 

0. 015 

0. 030 

1. 0 

1. 0 

0. 06 

0. 06 

2. 5 

6. 0 

2. 0 

2. 0 

0. 5 

0. 5 

0. 4 

0. 4 

0. 1 

0. 1 

0. 5 

0. 5 

2. 0 

2. 0 

2. 0 

2. 0 

0. 005 

0. 010 

0. 01 

0. 02 


M 

W, lb 

Center-of-gravity position, in. _ 

q, percent 

I y, percent 

a n , g units 

an 

Unswept-wing model, g units 

All other models, g units 

a, deg 

da n .nose 

> percent 

aCln,cg 

da n , 

percent 

P, sec 

T 1 / 2 , sec 


as and < ^ n ' 71086 > because these ratios determine the 
dot aa n% eg 

x 

accuracy of C L<x and respectively, rather than the abso- 
lute values of each quantity. The calculated accuracies of 
Cm q -\-C m6c are given in percent rather than as absolute mag- 
nitudes because the values of C mq -\-C ma have widely different 
numerical magnitudes for the three configurations flown in 
this investigation. 

Based on reference 13, the estimated possible error of each 
parameter is then the root-mean-square value of all the 
additive errors listed in table V, and these root-mean -square 
values are also given in table V. 

DISCUSSION 

SUMMARY OF RESULTS 

Lift. — The variations with Mach number of the total 
lift-curve slopes and the exposed-wing lift-curve slopes are 
shown in figure 27 for the four models as obtained from plots 
of C L against a as in figures 13 and 14. Data were obtained 
to lower Mach numbers for the unswept-wing model than 
for the other three models because this model did not pitch 
up at high angles of attack and execute violent lateral and 
longitudinal motions as did the other three models. 

The data for the delta-wing model (fig. 27 (a)) indicated 
nonlinearity with lift coefficient over the entire Mach num- 
ber range with the lift-curve slope of the total configuration 
at the high lifts being somewhat lower than that at the 
lower lifts. It is believed that this condition primarily 
resulted from the horizontal tail being in a position such 
that the downwash variation with angle of attack cl e/da 


TABLE V 

CALCULATED ACCURACY OF PARAMETERS 


[Increments may be positive or negative] 


Increments in 
parameters due 
to probable 
errors in the 
quantities 
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was greater as the lift increased. This type variation of 
del da with C L is indicated by the data in reference 14. The 
exposed-wing lift-curve slope (fig. 27 (a)) indicated only a 
small amount of nonlinearity over the lift range covered 
and supports the foregoing statement concerning the down- 
wash variation. Estimates of the exposed -wing lift-curve 
slope at supersonic speeds from the approximate linearized 
theory of reference 15 indicate very good agreement with 
the measured values. 

The lift-curve slopes for the two swept-wing models in 
figure 27 (b) illustrate the usefulness of the rocket-propelled- 
model test technique for the investigation of aeroelasticity. 
The variation of lift-curve slope with Mach number is 
shown in figure 27 (b) for the exposed wings and complete 
models. The slopes are taken near a lift coefficient of 0.30 
since data were available at this lift coefficient for both 
horizontal- tail deflections. The duralumin -wing lift-curve 
slopes for the complete model and exposed wing increase 
rapidly from a Mach number of 0.92 to 0.95 and then grad- 
ually decrease as the Mach number increases. The steel- 
wing lift-curve slopes show the same variation with Mach 
number but, as would be expected, have higher values than 
for the more flexible duralumin wing. Also shown in figure 
27 (b) is the lift-curve slope of the exposed rigid wing as 
obtained by the method described in the section entitled 
“Aeroelastic Analysis.” 

The lift data obtained from the unswept-wing configura- 
tion are given in figures 27 (c) and 28. The data in figure 
27 (c) indicate some nonlinearity of the lift curves in the 
region from M= 0.75 to M= 1.0 as evidenced by the different 
slopes obtained for lift coefficients in the regions C L = 0 and 
0.4. No values of C La could be determined near C L = 0 at 
supersonic speeds because of the small amplitudes of oscilla- 
tion following positive elevator deflections. The indications 
are that the nonlinearity disappears at Mach numbers above 
1.0. Figure 28 presents a summary of the lift information 
obtained on the unswept-wing model. At subsonic Mach 
numbers the model reached very high positive angles of 
attack each time the control was deflected in the negative 
direction. An examination of the telemeter record of normal 
acceleration indicated the onset of a high-frequency oscilla- 
tion of about 120 cycles per second, the wing first bending 
frequency. The value of lift coefficient where this oscillation 
began was taken to define the buffet boundary shown in 
figure 28. Beyond the buffet boundary the model apparently 
stalled. The points where this occurred are indicated by the 
maximum lift coefficient C Ltmax . These maximum lift coeffi- 
cients were obtained under dynamic conditions and are not 
necessarily the same as those obtained during static tests. 
The maximum rate of change of angle of attack in terms of 


the nondimensional factor 


da c 

dt 2V 


was 4.9 radians at a Mach 



(a) Delta-wing model. 

(b Swept-\\ ing models. Cl ^ 0.30. 
(c) Unswept-wing model. 

Figure 27. — Lift-curve slopes. 


number of 0.83 and the results of reference 16 indicate that 
any dynamic effect would be negligible at this test condition. 

Static stability. -The static stability characteristics of all 
models are given in figure 29 as the variation of aerodynamic- 
center position with Mach number. 

Results for the delta-wing model were obtained from both 


an analysis using the period of the oscillation and from two- 
accelerometer pitching-moment data; the agreement of the 
aerodynamic-center values determined from these two meth- 
ods is very good. The fairly smooth variation of static 
stability with Mach number was similar to variations ob- 
tained from the tailless delta-wing configurations reported 
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Figure 28. — Summary of buffet and maximum lift information. 
Unswept-wing model. 


in reference 1. The stability decreased with increase in lift 
coefficient. It is believed that the primary factor contribut- 
ing to the decrease in stability with increase in lift was an 
increase in del da at the tail with an increase in lift. This 
increase was associated with the high tail position and was 
mentioned in the preceding section of this report as a possible 
cause of the reduction in total-lift -curve slope with an increase 
in lift coefficient. 

The aerodynamic-center position for both swept-wing 
models is shown in figure 29 (b). The curves for both 
models were obtained by using the “period” method of 
analysis. The points for the steel-wing model were obtained 
from the double-differentiation method (eqs. (16) to (18)) 
and the points for the duralumin-wing model were obtained 
from the two-accelerometer method. With the exception of 
the one point at M= 1.27 for the steel-wing model, the 
agreement between the points and the curve is good. 

The aerodynamic-center position for the steel-wing model 
has the same variation with Mach number as the duralumin- 
wing model over the Mach number range 1.02 and 1.27. 
Over this Mach number range, however, the aerodynamic- 
center position is 3 to 6 percent more rearward than that 
for the steel-wing model. From the strip-theory method 
of the section entitled “Aeroelastic Analysis” the forward 
movement in aerodynamic-center position due to the in- 
board movement of the span load was calculated for the 
duralumin-wing model. A comparison between the calcu- 
lation made for a 0.50-chord loading and the measured 
difference in aerodynamic-center position for the steel- and 
duralumin-wing models in figure 29 (b) shows that nearly all 
the loss in stability of the duralumin -wing model may be 
accounted for by the inboard movement of load on the wing. 

The aerodynamic-center location for the unswept-wing 
configuration (fig. 29 (c)) indicates some nonlinearity in the 
pitching-moment curves at subsonic speeds. The curve 
faired through the data for C L =0A at subsonic speeds has 
been dotted to indicate that the results were obtained from 
measuring periods during oscillations in which the angle-of- 
attack vane was against a stop during part of the oscillation. 



(a) Delta wing model, 
(b) Swept-wing models, 
(c) Unswept-wing model. 


Figure 29.— Aerodynamic-center locations. 


As the Mach number is increased above about 0.82, the 
aerodynamic center first moves forward to its most forward 
location at 35 percent of the mean aerodynamic chord at a 
Mach number of 0.90. It then moves rearward at higher 
Mach numbers. 

Damping in pitch. — The variation with Mach number of 
the damping-in-pitch parameter Cm q + C m6t for all four models 
is shown in figures 22 and 30. These results were obtained 
from application of equation (31) by using the lift-curve-slope 
data of figure 27 along with the time to damp to one-half 
amplitude as discussed and illustrated previously. 
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For the delta-wing configuration an increase in the nega- 
tive value of Cm q +Cm a with increase in lift coefficient is in- 
dicated. (See fig. 22.) An increase in de/da with increase in 
lift coefficient mentioned previously could account for this 
increase in C mq +C ma through its effect in increasing C m ^ 

The values of C mq +C ma for both swept-wing models are 
shown in figure 30 (a). The duralumin-wing-model data are 
for a eenter-of-gravity position 20 percent of the mean aero- 
dynamic chord more forward than that for the steel-wing 
model. This difference in center-of-gravity position was cal- 
culated to have a small effect on the damping-in -pitch deriva- 
tive. The damping in pitch derivative for the duralumin- 
wing model shows a rapid variation with Mach number in the 
transonic speed range. This same variation might have 
occurred for the steel-wing model also, but damping data 
were not obtained for the steel-wing model below a Mach 
number of 1.0. The lower value of C mq -\-C m6t for the dural- 
umin-wing model from Mach numbers 1.02 and 1.2 may not 
be due to aeroelasticity and is more likely to be an indication 
of the accuracy of the damping-in-pitch derivative. 

The results for the unswept-wing configuration are given in 
figure 30 (b). These results indicate a large increase in the 
negative value of C mq +C ma in the transonic range as was the 
case for the duralumin swept-wing model; however, the 
variation with Mach number after the peak is not as abrupt. 

These results presented for C mq + C ma indicate the wide 
range of numerical values associated with this parameter; 
this wide range makes it impossible to assess the dynamic 
stability of a configuration by means of this parameter alone. 

Drag. — There are several ways of presenting the drag data 
obtained. The drag polars for the variation of C D with C L 



(a) Swept-wing models. 

(b) Unswept-wing model. 

Figure 30. — Pitching-moment damping factor. 


were illustrated for the swept-wing configuration. (See fig. 
15.) Figures 31 to 33 give the drag characteristics of all four 
models as a function of Mach number and illustrate some 
methods of plotting the results. 

The drag data for the delta-wing configuration is given as 
the variation of drag coefficient with Mach number at several 
constant lift coefficients (fig. 31) and at two control deflec- 
tions. Plots such as this give a rather complete picture of the 
drag characteristics of a configuration. For instance, the 
drag increment due to control deflection at lift coefficients of 
0.2 and 0.3 is very evident. A large portion of this change in 
drag coefficient resulting from a change in tail incidence is the 
change in the streamwise component of the tail normal-force 
coefficient with change in tail incidence. 

Another way of presenting drag data is illustrated by the 
results for the swept-wing configuration. The drag data for 
the steel- and duralumin-wing models of the swept-wing con- 
figuration are given in figure 32 as the variation of C Dtm in ail( l 


dC D 


with Mach number. The minimum drag values for the 


horizontal-tail setting of 0.14° for the steel-wing model and 
— 0.72° for the duralumin-wing model are in excellent agree- 
ment, the accuracy of the measurements and the small 
difference in horizontal-tail deflection being considered. 



Figure 31. — Variation of drag coefficient with Mach number and lift 
coefficient for the delta-wing model. 
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The drag due to lift of the models is presented in figure 
32 (b) for lift coefficients of 0.50 and in figure 32 (c) for lift 
coefficients of 0.15 since the induced drag was found to be 
nonlinear with lift coefficient. Also shown in figures 32 (b) 
and 32 (c) are the ideal induced-drag factor \/ttA and the 
drag due to lift for no leading-edge suction 1/57.3 C L<x . For 
the 1/57.3 C La comparison, the lift-curve slopes were meas- 
ured near the trim lift coefficient corresponding to elevator 
deflections of —3.0° and 0.14° for the steel-wing model and 
— 4.60° and —0.72° for the duralumin-wing model. The 
lower induced drag for the duralumin-wing model suggests 
the possibility of a more efficient span loading for the flexible 
wing as a result of the inboard movement of the load when 
the wing deflects. 

The results for the unswept-wing configuration are pre- 
sented in figure 33. The longitudinal accelerometer in this 



M 


(a) Minimum drag. 

(b) Drag due to lift Cl « 0.50. 

(c) Drag due to lift Cl « 0.15. 

Figure 32. — Drag characteristics of the swept-wing models. 


model was calibrated to cover a sufficient range to include 
the accelerations during booster burning (about 1 8g) and 
thus did not give very good accuracy on the accelerations 
(about —\g to — 4gf) developed during the time the drag data 
were obtained. This is evidenced by the scatter of data in 
figure 33 (a). The minimum drag at supersonic speeds (fig. 
33 (a)) is fairly high and is due mostty to the fuselage which 
is not a particularly good shape for supersonic speeds. 

The effect of lift on drag is shown in various ways in figures 
33 (b) and 33 (c). The (. L/D) ma£ (fig. 33 (b)) decreases by 
about one-half as the Mach number increases from 0.8 to 
1.0. This decrease is apparently due to the increase in mini- 


mum drag because 


dC D 

dC L 2 


does not increase in this Mach 



(a) Minimum drag. 

(b) Maximum lift-drag ratio. 

(c) Drag due to lift. 

Figure 33. — Drag characteristics of the unswept-wing model. 
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number region, as shown in figure 33 (c). The dashed curve 


in figure 33 (c) is a plot of 


n which should equal the 
57 L(x 


value of ( \[ d - 2 for a wing with the resultant force normal to 

the chord plane. The agreement between this curve and 

the measured values of is good except at the highest 

a( L 


Mach numbers. 

Control effectiveness. -The longitudinal control con- 
sisted of all -movable stabilizers which were the same on all 
models except that on two models the gaps at the roots 
were sealed as mentioned previously. When based on the 
same dimensions, the values of C ms should thus be the same 
for all models except for the effect of the gaps. Figure 34 
shows values of C ms for the three configurations multiplied 
by Sc /l t in each case to furnish a common basis for compari- 
son. The single curve shown for each configuration is the 
average of the values obtained by the methods illustrated 
previously. (See fig. 23 (a).) The cross-hatched area indi- 
cates the region which contained all the individual values 
of C ms as determined by all the methods. 

When the nonlinearity of the pitching-moment data dis- 
cussed earlier is considered, the agreement of the data shown 
in figure 34 is believed to be reasonable since only two con- 
trol deflections were available for each model and an assump- 
tion of linearity was thus required. Any effect of the gap 
at the horizontal- tail root on the control effectiveness is 
masked by the inaccuracy of determining C ms for the models 
of this investigation. 


In the speed region under consideration, the Mach num- 
ber is, of course, the most important similarity parameter 
and is used as the basic comparison variable. Viscous- 
flow effects are dependent upon Reynolds number and 
figure 35 shows the Reynolds numbers for the free-flight 
models and for two typical full-scale airplane flight condi- 
tions. The model Reynolds numbers approach those for 
airplanes at altitudes that are operationally feasible for 
these configurations and Mach numbers. 

Experimental and theoretical analyses of unsteady aero- 
dynamic effects indicate that these effects are dependent 

on the nondimensional frequency parameter ~~ For the 

models discussed herein, the range of values of was from 

0.008 to 0.018. These values are considerably below the 
frequencies at which unsteady flow effects generally became 
appreciable and, in addition, are reasonably close to the 
probable values for full-scale airplanes of similar con- 
figuration. 

The dynamics of the motion are affected by several 
similitude factors. If equation (33) for C ma is rewritten in 
terms of nondimensional quantities, the usually small 
damping term being neglected, the following expression for 
the nondimensional frequency parameter is obtained: 



APPLICATION TO FULL-SCALE AIRPLANES 

When the aerodynamic characteristics have been de- 
termined from the model flight analysis, the performance 
and flying qualities of a similar full-scale airplane may, of 
course, be determined by substituting these quantities 
along with the airplane mass and geometric characteristics 
into the equations of motion. Certain questions of scaling 
arise when this procedure is used and some consideration 
has been given to this subject. 



Figure 34. — Comparison of control-effectiveness data. Average 
values used to determine curves. 



Figure 35. — Reynolds numbers based on mean aerodynamic chord for 
models and typical values for an airplane with c= 10 ft. Condition 
1, 60,000 feet. Condition 2, 40,000 feet. 
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Values of the relative density factor /z for the models of the 
present investigation are shown in figure 36. Figure 36 
also has an auxiliary scale that shows the equivalent altitude 
for an airplane with a chord of 10 feet and a wing loading of 
100 pounds per square foot. It appears that the conditions 
under which the models were flown approximate those for 
typical airplanes at reasonable altitudes. A comparison 
of the values of k Y /c for the present models and several 
current full-scale airplanes indicates that the values for the 
models are appreciably larger. This effect is at least partially 
compensated for in most cases by further forward locations 
of the center of gravity on the models, resulting in larger 
values of C ma . If dynamic similarity of the oscillation fre- 
quency only is desired, several factors are thus available for 
adjustment during test planning such as altitude, weight, 
radius of gyration, and center-of-gravity position. 

The conditions governing the damping are more restrictive 
than those applying to the frequency. By making use of a 
nondimensional damping factor, such as the number of 
cycles required to damp to one-half amplitude, and deriving 
this factor from equations (26), (27), (29), and (50), the 
following expression can be obtained: 


r - • 312 V~ <7 ’”^ 

)2 C lc Y Cmq+C^ 


(51) 


Equation (51) shows that k Y /c must be the same for the model 
and the full-scale airplane in order for the proportions of 
the damping contributed by the C L<X and C m to be the 

same. The larger values of k Y /c generally occurring on the 
models, as mentioned previously, results in a larger contri- 
bution of the C La - term to the damping as compared with 
typical full-scale conditions. This condition has a dis- 



Figure 36. — Relative density factors for models and equivalent values 
of altitude for an airplane with c=10 ft and JU/N^lOO lb/sq ft. 


advantage, in one sense, because the accuracy of determin- 
ing Cm q +Cm a is correspondingly decreased in the model 
tests. An important advantage from this increased effect 
of C L<x may exist in some cases, however, because damped 
model motions may be obtained to permit determination of 
unstable values of which could lead to dynamic 

instability under full-scale conditions. 

In general terms, it can be stated that the motions ob- 
tained on the models can be taken as applying to a full-scale 
airplane having the same center-of-gravity position, relative 
density, and nondimensional radius of gyration. The time 
histories of such motions may be applied directly to full-scale 
conditions if the time scale of the motions is multiplied by the 
ratio of the airplane size to the model size and the ratio of 
speed of sound at model altitude to speed of sound at airplane 
altitude and if the scale of angular velocities is reduced by 
the same factor. The latter statement applies to nonlinear 
motions as well as linear motions and also applies to com- 
bined pitching, rolling, and yawing motions if all the non- 
dimensional radii of gyration are the same on the model and 
the full-scale airplane. 

A further factor that needs to be considered is the scaling 
of aeroelastic effects. The parameter (iS(0lL) ref utilized in 
figure 25 to relate data from different stiffness wings is one 
scaling factor that may be used. The range of the static- 
pressure ratio obtained on all models of this report is given 
in figure 37 and values of (i 9/L ) for the duralumin swept 
wing are given in figure 10. Because of the low altitude and 
resulting high dynamic pressures of the model test, the 
aeroelastic effects observed on the solid duralumin wing are 
directly applicable to an airplane with a more flexible wing 
at a higher altitude. 

In case the model results are not directly applicable to the 
full-scale airplane, it is possible to correct the model results 
to rigid conditions, as illustrated herein, and then to use 
these results with the airplane elastic characteristics to 
determine the effects of flexibility on the airplane. 

CONCLUDING REMARKS 

A method has been set forth in some detail for determining 
aerodynamic performance and longitudinal stability param- 
eters from free-flight tests of rocket-propelled models. The 
feasibility of the method has been demonstrated by data 
measured through the transonic speed range on several air- 
plane configurations. The data are obtained from analysis 
of the model motions following control-surface disturbances. 
Data can be obtained over a considerable Mach number 
range during each flight. 



28 


REPORT 1337— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


The method makes use of accelerometers and an angle-of- 
attack indicator to permit instantaneous measurements of 
lift, drag, and pitching moments. Thus, nonlinear charac- 
teristics can be studied. For configurations which do not 
exhibit excessive lateral motions at high angles of attack, 
the maximum lift and buffeting characteristics may be deter- 
mined. Effects of aeroelasticity may be readily obtained. 
Additional information can be obtained by analysis of the 
characteristics of the short-period free oscillations resulting 
from rapid control-surface deflections, but this procedure re- 
quires the assumption of linearity of the aerodynamic quan- 
tities during each oscillation. Some measure of nonlineari- 
ties may be obtained with the latter method, however, by 
studying oscillations at different trim angles of attack. Al- 
though the models decelerated continuously during the data- 
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